A Higher Order Scheme for a Tangentially Stabilized Plane Curve Shortening Flow with a Driving Force

We introduce a new higher order scheme for computing a tangentially stabilized curve shortening flow with a driving force represented by an intrinsic partial differential equation for an evolving curve position vector. Our new scheme is a combination of the explicit forward Euler and the fully-implicit backward Euler schemes. At any discrete time step, the solution is found efficiently using a few semi-implicit iterations. Basic properties of the new scheme are proved in the paper and its precision is tested by comparing the results with known analytical solutions. For any choice of the time step, the new higher order scheme gives exact radius of evolving uniformly discretized circles in case of flow by curvature and in case of rotation by a constant tangential velocity. Such properties do not hold for other schemes solving flow by mean curvature like the classical explicit, semi-implicit or fully-implicit schemes. In general, the scheme is second order accurate, which is shown by comparing a numerically evolving encompassed area with know analytical expression. The behavior of the scheme is discussed on representative examples and its advantages with respect to the balance between CPU time and precision are shown.